f[x_] := (1/(Exp[((x - μ)/(kB*T))] + 1)); h = -D[f[x], T]; 
    rule1 = {kB -> 1, μ -> 1};
f1 = f[x] /. rule1;
p1 = 
      Plot[Evaluate[
          Table[f1, {T, 0.002, 0.02, 0.002}]], 
{x, 0.9, 
          0.99}]

f[x_] := (1/(Exp[((x - μ)/(kB*T))] + 1)); h = -D[f[x], T]; 
    rule1 = {kB -> 1, μ -> 1};
f1 = f[x] /. rule1;
p1 = 
      Plot[Evaluate[
          Table[f1, {T, 0.002, 0.02, 0.002}]], 
{x, 1.01, 
          1.10}]

f[x_] := (1/(Exp[((x - 1)/T)] + 1));
\(p1 = 
      Plot[Evaluate[
          Table[f[x], {T, 0.002, 0.02, 0.002}]], 
{x, 0.9, 
          1.10}, PlotRange -> {{0.9, 1.10}, {0, 
              1}}, 
PlotStyle -> 
          Table[Hue[0.1  i], {i, 1, 10}], 
Prolog -> 
          AbsoluteThickness[2], 
        AxesLabel -> {"E", "F(E)"}, 
Background -> GrayLevel[0.7]]\)

g1 = D[f[x], x]; 
    p2 = Plot[
        Evaluate[Table[-g1, {T, 0.002, 0.02, 0.002}]], {x, 0.9, 1.10}, 
        PlotRange -> {{0.9, 1.10}, {0, 130}}, 
        PlotStyle -> Table[Hue[0.1  i], {i, 1, 10}], 
        Prolog -> AbsoluteThickness[2], 
        Background -> GrayLevel[0.7]]

g[
          x_] := ((Sqrt[x])/(Exp[((x - 1)/(1  T))] + 1));
\(p1 = 
      Plot[Evaluate[
          Table[g[x], {T, 0.001, 0.05, 0.01}]], 
{x, 0, 
          1.30}, PlotRange -> {{0, 1.5}, {0, 
              1}}, 
PlotStyle -> 
          Table[Hue[0.2  i], {i, 1, 10}], 
Prolog -> 
          AbsoluteThickness[2], 
        AxesLabel -> {"E", "D(E)f(E)"}, 
Background -> GrayLevel[0.7]]\)

f[x_] := (1/(Exp[((x - 1)/T)] + 1));
\(p1 = 
      Plot[Evaluate[
          Table[f[x], {T, 0.002, 0.02, 0.002}]], 
{x, 0.9, 
          1.10}, PlotRange -> {{0.9, 1.10}, {0, 
              1}}, 
PlotStyle -> 
          Table[Hue[0.1  i], {i, 1, 10}], 
Prolog -> 
          AbsoluteThickness[2], 
        AxesLabel -> {"E", "F(E)"}, 
Background -> 
          GrayLevel[0.7]]\)
ImageSize -> 0.5

g[
          x_] := ((Sqrt[x])/(Exp[((x - 1)/(1  T))] + 1));
\(p1 = 
      Plot[Evaluate[
          Table[g[x], {T, 0.001, 0.05, 0.01}]], 
{x, 0, 
          1.30}, PlotRange -> {{0, 1.5}, {0, 
              1}}, 
PlotStyle -> 
          Table[Hue[0.2  i], {i, 1, 10}], 
Prolog -> 
          AbsoluteThickness[2], 
        AxesLabel -> {"E", "D(E)f(E)"}, 
Background -> 
          GrayLevel[0.7]]\)
ImageSize -> 0.5

f[x_] := (1/(Exp[((x - 1)/T)] + 1));
\(p1 = 
      Plot[Evaluate[
          Table[f[x], {T, 0.002, 0.02, 0.002}]], 
{x, 0.9, 
          1.10}, PlotRange -> {{0.9, 1.10}, {0, 
              1}}, 
PlotStyle -> 
          Table[Hue[0.1  i], {i, 1, 10}], 
Prolog -> 
          AbsoluteThickness[2], 
        AxesLabel -> {"E", "f(E)"}, 
Background -> 
          GrayLevel[0.7]]\)
ImageSize -> 0.5

ImageSize -> 0.5`

g[
          x_] := ((Sqrt[x])/(Exp[((x - 1)/(1  T))] + 1));
\(p1 = 
      Plot[Evaluate[
          Table[g[x], {T, 0.001, 0.05, 0.01}]], 
{x, 0, 
          1.30}, PlotRange -> {{0, 1.5}, {0, 
              1}}, 
PlotStyle -> 
          Table[Hue[0.2  i], {i, 1, 10}], 
Prolog -> 
          AbsoluteThickness[2], 
        AxesLabel -> {"E", "g(E)f(E)"}, 
Background -> 
          GrayLevel[0.7]]\)
ImageSize -> 0.5

f[x_] := (1/(Exp[((x - 1)/T)] - 1));
\(p1 = 
      Plot[Evaluate[
          Table[f[x], {T, 0.002, 0.02, 0.002}]], 
{x, 0.9, 
          1.10}, PlotRange -> {{0.9, 1.10}, {0, 
              1}}, 
PlotStyle -> 
          Table[Hue[0.1  i], {i, 1, 10}], 
Prolog -> 
          AbsoluteThickness[2], 
        AxesLabel -> {"E", "F(E)"}, 
Background -> GrayLevel[0.7]]\)

f[x_] := (1/Exp)[((x - 1)/T)];
\(p1 = 
      Plot[Evaluate[
          Table[f[x], {T, 0.002, 0.02, 0.002}]], 
{x, 0.9, 
          1.10}, PlotRange -> {{0.9, 1.10}, {0, 
              1}}, 
PlotStyle -> 
          Table[Hue[0.1  i], {i, 1, 10}], 
Prolog -> 
          AbsoluteThickness[2], 
        AxesLabel -> {"E", "F(E)"}, 
Background -> GrayLevel[0.7]]\)

g[
          x_] := ((Sqrt[x])/(Exp[((x - 1)/(1  T))] + 1));
\(p1 = 
      Plot[Evaluate[
          Table[g[x], {T, 0.001, 0.05, 0.01}]], 
{x, 0, 
          1.30}, PlotRange -> {{0, 1.5}, {0, 
              1}}, 
PlotStyle -> 
          Table[Hue[0.2  i], {i, 1, 10}], 
Prolog -> 
          AbsoluteThickness[2], 
        AxesLabel -> {"E", "N(E)=D(E)f(E)"}, 
Background -> 
          GrayLevel[0.7]]\)
ImageSize -> 0.5

g[
          x_] := ((Sqrt[x])/(Exp[((x - 1)/(1  T))] + 1));
\(p1 = 
      Plot[Evaluate[
          Table[g[x], {T, 0.001, 0.05, 0.01}]], 
{x, 0, 
          1.30}, PlotRange -> {{0, 1.5}, {0, 
              1}}, 
PlotStyle -> 
          Table[Hue[0.2  i], {i, 1, 10}], 
Prolog -> 
          AbsoluteThickness[2], 
        AxesLabel -> {"E", "N(E)=g(E)F(E)"}, 
Background -> 
          GrayLevel[0.7]]\)
ImageSize -> 0.5

g[
          x_] := ((Sqrt[x])/(Exp[((x - 1)/(1  T))] + 1));
\(p1 = 
      Plot[Evaluate[
          Table[g[x], {T, 0.001, 0.05, 0.01}]], 
{x, 0, 
          1.30}, PlotRange -> {{0, 1.5}, {0, 
              1}}, 
PlotStyle -> 
          Table[Hue[0.2  i], {i, 1, 10}], 
Prolog -> 
          AbsoluteThickness[2], 
        AxesLabel -> {"E", "N(E)=g(E)f(E)"}, 
Background -> 
          GrayLevel[0.7]]\)
ImageSize -> 0.5